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Abstract. - We examine the non-extensive approach to the statistical mechanics of Hamiltonian 
systems with H — T + V where T is the classical kinetic energy. Our analysis starts from the basics 
of the formalism by applying the standard variational method for maximizing the entropy subject 
to the average energy and normalization constraints. The analytical results show (i) that the non- 
extensive thermodynamics formalism should be called into question to explain experimental results 
described by extended exponential distributions exhibiting long tails, i.e. g-exponentials with q > 1, 
and (ii) that in the thermodynamic limit the theory is only consistent in the range < q < 1 where 
the distribution has finite support, thus implying that configurations with e.g. energy above some 
limit have zero probability, which is at variance with the physics of systems in contact with a 
heat reservoir. We also discuss the ((/-dependent) thermodynamic temperature and the generalized 
specific heat. 



Introduction. — The observation of natural phenomena and of laboratory experiments 
provides a wide spectrum of experimental results showing distributions of data that deviate 
from exponential decay as predicted for Boltzmann behavior [T]. It was the goal of non- 
extensive statistical mechanics developed originally by Tsallis [3] to offer a new approach 
to explain the non-Boltzmann behavior of non-equilibrium systems 0- More precisely the 
primary motivation for non-extensive thermodynamics is as a way to understand deformed 
exponential distributions (such as g-exponentials exhibiting long tails when q > 1) found 
empirically in many areas of physics and other scientific disciplines [S] . The interest raised by 
this new approach has grown over the years and has produced an abundant literature [5] re- 
flecting new theoretical developments and a considerable number of applications to subjects 
as diverse as defect turbulence, energy distribution in cosmic rays, earthquake magnitude 
value distributions and velocity distributions in micro-organism populations or distributions 
of financial market data. Parallel to these developments, critical analyses were presented 
questioning the merits of the non-extensive theory [7] and thereby of its applications as well. 
These criticisms raised questions that often gave rise to ontological conflicts [8] . Most of the 
criticisms are based on phenomenological analyses and thermodynamic arguments question- 
ing the compatibility of the theory with classical statistical thermodynamics. More recently, 
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1 We note that other approaches exist such as super-statistics [3] and the fractional and nonlinear Fokker- 
Planck equations (see e.g. [4] and references therein) 
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the non-extensive theory was also re-examined on the basis of analyses demonstrating the 
necessary discreteness of systems to which the theory applies [9] and the limits of validity 
of the non-extensive formalism for a Hamiltonian system, the g-ideal gas, a model system of 
independent quasi-particles where the interactions are incorporated in the particles proper- 
ties [TU1[TT] . Here we merely adopt the viewpoint of analytical rigor to establish the validity 
limits of the non-extensive formalism for the general class of classical Hamiltonian systems 
with continuous variables and consequently of the class of physical systems to which the 
non-extensive interpretation applies. 

In the formalism, the q-exponential distributions arise as the result of a variational 
calculation maximizing the generalized entropy, the g-entropy, under the constraints of nor- 
malizability of the distribution function and of a prescribed average energy [5]. The goal 
of the present work is to investigate the exact form of the distribution so derived for clas- 
sical Hamiltonian systems from both the usual Tsallis entropy [5] and the homogeneous 
entropy [T3]. The conclusion is found to be the same in both cases: the theory is only 
consistent in the thermodynamic limit for < q < 1. For finite systems of N particles, the 
upper limit is 1 + a/N for a ~ 0(1). This means (i) that the non-extensive thermodynamics 
formalism cannot be used, at least in any straightforward way, to explain phenomena for 
which one observes that q takes a value q > 1 (corresponding to asymptotically power-law 
distributions) i.e. when the distributions exhibit extended exponential forms with long tails, 
and (ii) that, within the validity domain, the distribution has finite support, thus imply- 
ing that configurations with e.g. energy above some limit have zero probability, a strange 
situation for systems with typical molecular potentials which are steeply repulsive at small 
distances. We treat successively the case according to the development based on the Tsallis 
entropy S q and the case based on the homogeneous entropy . The reason for examining 
the two cases is the criterion of stability against perturbations of the probability distribu- 
tion function, or Lesche stability |12j : the homogeneous entropy was shown to be Lesche 
stable while the Tsallis entropy is not [T3]. We also present the results for the generalized 
(g-dependent) thermodynamic properties of Hamiltonian systems in both cases. 

Tsallis entropy. — Non-extensive statistical mechanics is developed on the basis of 
three axioms: (i) the g-entropy for systems with continuous variables is given by [B] 

s.-*. 1 - (1) 

where T is the phase space variable and K must be a quantity with the dimensions of 
[r] 9_1 , i.e. K = fi 01 ^^ 1 } (£) denotes the dimension of the system and TV the number of 
particles) and the classical Boltzmann-Gibbs entropy is retrieved in the limit q — > 1; (ii) the 
distribution function p (r) is slaved to the normalization condition 

1 = J P (T)dV- (2) 

(iii) the internal energy is measured as 

f f P q (T)HdT 

U = J P ^ HdT = fpl(r)aT ' < 3 > 

where P q (T) is the escort probability distribution [14] which is the actual probability mea- 
sure. We consider Hamiltonian systems with H = T + V where T is the classical kinetic 
energy. The distribution function p (T) is obtained by maximizing the g-entropy subject to 
the normalization © and average energy ([3]). Introducing the Lagrange multipliers a and 
7, the variational method leads to 
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which is solved to yield 



p(T) = (l-(l-g)7(H-EQ)^ ^ (5) 



J (l-(l-q)j(H-U))^dT 



with 7 = k j p q(Y)dr anc ^ where the normalization condition ([3]) has been used to eliminate 
the multiplier a. The other, 7, is determined from the energy constraint (j3J); using in 
fl3J), we have 

I (1 - (1 -q) 7 (H- U))^ (H- U)dY 

= I 7 , (6) 

(1 - (1 - q) 7 (H - U))^ dT 



which by multiplying the numerator by (1 — (7)7 and adding and subtracting one gives 

(1 - (1 - q) 7 (H - U))^ dT= [ (1 - (1 - q) 7 (H - IT)) 1 * 5 dT . (7) 



Note that the sign of the factor f(T) = 1 — (1 — q) 7 (H — U) is, so far, arbitrary so that 
we allow for the cancellation of a (possible imaginary) factor throughout these equations. 
However, this is only possible in the expression for the distribution <j5j) if the sign of f(T) 
is independent of T. Given this fact, the sign can be fixed by making the substitution 

q 1 

/(F) = s|/(r)| in the energy constraint which requires that s 1 -' = s 1 -? so that s = 1. 
Since the kinetic energy is in principle unbounded, the constant sign of the argument means 
that the distribution may have to be restricted to some subset of phase space so that it 
should be written as 

(r , (I - (I - q)j(H - U))^ 0(1 - (1 - q) 7 (H - U)) 

P\ L ) = —f ; ) (°) 



(l-(l-q) 1 (H-U))—G(l-(l-q) 1 (H-U))dT 

where 9(.t) is the step function which is one for x > and zero otherwise. We note that the 
introduction of the step function is in fact a redefinition of the variational problem in the 
sense that we have replaced p(T) with p(T)Q(f(T)) in Eq.([T][3j) and then maximized with 
respect to p(T). If we do not do this, then there is simply no solution to the variational 
problem which is real and no n- negative for all T. To see this, we note from the variational 
equation (01 that if p(T) vanishes for some value of T then it necessarily follows that a — 0, 
which is untenable since a will generally assume a non-zero value due to the normalization 
condition. Therefore, we can only restrict the support of the distribution by redefining the 
variation problem. We now turn to the determination of the possible values for q. 

Case: q < 1. The exponent occuring in the distribution is > so that if 7 > 0, 
then the distribution must have finite support since for some sufficiently large value of T , 
the argument, f(T), becomes zero, and from (|5|) we have 

cxp g (-7(ff-EQ) 
M j ~/exp 9 (- 7 (H-t/))dr' W 

where exp q (x) = (1 + (1 — g)i) Tr '6(l + (1 — q)x) is the g-exponential function. 

If 7 < 0, then /(r) is always positive, and the distribution function never goes to zero 
so that there can be no finite support. This however leads to another problem since f(T) is 
unbounded as the kinetic energy increases which, in turn, means that the integral of f(T) 
over momenta will diverge so that the distribution cannot be normalized. We conclude that 
7 < is not allowed. We remark that it might be thought that one could introduce a limit 
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on the kinetic energy, but this is not in keeping with the proposal that the nonextensive 
distribution be a generalization of the canonical distribution which describes open systems 
in contact with a reservoir. 

Case: q > 1. The exponent being then < 0, we write the distribution function as 

(r) _ (i + |i-g|7(g-t/))' lT ^ | e(i + |i- g | 7 (g-t/)) (10) 

J (1 + |1 - g|7(JT - U))~^ 6 (1 + |1 - q\ 7 (H - U)) dT 

If 7 > 0, the distribution will be normalizable provided that (J^J 5 ?) l q ^ IS integrable over 

1 I 

d ND p (for large p) which is to say that (P 2 ) 1 ~ q P ND ~ 1 dP be integrable for P — > 00; this 

means we must have — 1 > ND — 1 — 2 , or 1 < q < 1 + which condition reduces 

to the classical Boltzmann result (q = 1) in the thermodynamic limit 0. 
The other possibility is 7 < 0. We then write the distribution as 

(r) = (1 - M |i - gj Cg - 1/))" 1 ^ 1 e (1 - |tJ |i - q \ (h - u)) 

\l - | 7 | |1 - q\ (H - U))~\t^\ 9 (1 - | 7 | |1 - q\ (H - U)) dT 

When integrated over momenta, this expression would show a singularity at T = r^ny^ + 
U — V, unless 

j(x- j: lP 2 y^d P m ~ (x - p^\^\ P 3ND-i dP „ j-x {x _ Y) \^- q \ Y mj^i dY 

is finite; this requires that 1 — 



l-q 



= 1 - > 0, i.e. q > 2 . But p q (T) must also 

be integrable and in order that the singularity be integrable imposes 1 — > which is 
incompatible with the condition q > 2. 

Furthermore the maximization condition demands that the second derivative of the g- 
entropy (JlJ be — ks K q p q ~ 2 < 0, which is satisfied if q > 0. So in summary, the distribution 
function exists in the thermodynamic limit for 7 > and < q < 1, and p (T) has the form 
of a g-exponential with finite support. 

Homogeneous entropy. — The homogeneous entropy was proposed as an alternative 
to the Tsallis entropy because it has various desirable properties that the Tsallis entropy 
does not share such as being Lesche stable and giving positive-definite heat capacities |13j . 
It is therefore interesting to ask whether it is subject to the same limitations as found for 
the Tsallis entropy. The analysis follows essentially the same lines as above except that in 
this case p (T) is the physical probability [T3] and so the energy is computed with the normal 
average. The homogeneous entropy reads |13j 

i-(Kf P HT)drY 

Sf = k B ^- >— , 12 

l-q 

where K is a quantity with the dimensions [r] q , i.e. K = h 1 , and the normalization 
and energy constraints are 



1 = J p{T)dT ; U = jp{T)HdT. 



Following a similar analysis as given in the previous section, the conclusions are that the 
condition for normalizability is < q < 1 and 7 > in which case the distribution function 
reads 



2 Note that p q (T) must also be integrable which imposes that —1 > ND — 1 — or q < 1 
but this condition is weaker than the constraint 1 < q < 1 + -^jr . 
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P(r) 



(exp g (- 7 (g -£/)))« 

far (cxp 9 (- 7 (H-u))Y 



(13) 



- where 7 is the Lagrange multiplier used to fix the average 



Here 7 = -, j— 

energy and p (T) has finite support. 

Thermodynamic properties. — Given that the formalism is constructed solely on 
the basis of the three axioms ([1]), ^ and ([3]), the consistent way to define the thermodynamic 
temperature is through the thermodynamic definition dS/dU = l/T. Considering q < 1 
and 7 > 0, we obtain from the Tsallis entropy (H} 



t t = 
g 



1 



fc B 7 



R— / exp (-7 (H - U)) dT 



q-l 



and from the homogeneous entropy (|12[) 

<1 



-iH 



cxp„ 



(H — U) \ dT 



l-q 



(14) 



(15) 



With the explicit expressions of 7 (see © and (flU)) ). (|15[) gives 7 = l/(k B T^), the analog 
of the classical expression, whereas the equivalent relation for the Tsallis temperature is 
only obtained in the limit q — > 1. 



Correspondingly, the expressions for the specific heat Cy — ( Jp- ) arc given by 



C 



ifBLY 1 

1 V T / A' 2 / (9 2 '!- 1 (r)(/3 T (ff 



" 2 C 1 - 



(16) 



with the classical notation /3 T = l/(fcgT^ ), and by 



9 9 



7 



fee 
~K 



p^r (r) (/3 ff (F-C/)) 2 (ir 



where P H = l/(k B T"), or 



C 



H 



k B (T*)* 



J dT P (T) {H-Uf C q (T) 



1 p 1 (r) 



(17) 



(18) 



with Cq(r) = qi \ — k ■ * s ^ ne generalization of the classical expression of 

the specific heat given in terms of the energy fluctuations: Cy = ((AE) 2 ) / (k B T 2 ). The 
thermodynamic temperatures are both positive while the generalized specific heat is only 
positive-definite in the case of the homogeneous entropyH 

Concluding comments. — Non-exponential distributions are widely observed in na- 
ture. Non-extensive thermodynamics was motivated, in part, as a means of explaining the 
origin of such distributions which arise naturally as a result of maximizing the generalized 
entropy with the usual constraints of the normalization of the distribution and of fixed 
average energy. We have shown that when this procedure is applied to Hamiltonian sys- 
tems, the resulting distributions only exist for the restricted range of < q < 1 + O (4>) . 



3 It was indeed shown that in the Tsallis formulation the specific heat can be negative 1101 . 
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Since the so-called "fat-tailed" distributions correspond to g > 1, this means that gen- 
eralized thermodynamics cannot be seen as an explanation of their occurrence for these 
systems. The problem with larger values of q has to do with the existence of the integrals 
over momenta due to the unboundedness of the kinetic energy. One way around this would 
be to redefine the formalism so as to restrict the momenta a priori by making the ansatz 
p (r) = (T* — T)p*(r), for some fixed positive number throughout the variational 
problem and maximizing to determine p*(T). However, this is obviously quite artificial and 
ad hoc since, for example, one could replace the step function by any function of T that goes 
to zero sufficiently quickly as T grows. This suggests the more straightforward conclusion 
that in the case of classical Hamiltonian systems, nonextensive thermodynamics does not 
provide a simple, natural explanation of distributions with fat tails. 

* * * 

This work was partly supported by the European Space Agency under contract number 
ESA AO-2004-070. 
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